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ABSTRACT 

A previous work (Monteserm et aUl2QQ8^ estimated the CMB variance from the three- 
year WMAP data, finding a lower value than expected from Gaussian simulations 
using the WMAP best-fit cosmological model. We repeat the analysis on the five-year 
WMAP data using a new estimator with lower bias and variance. Our results confirm 
this anomaly at higher significance, namely with a p- value of 0.31%. We perform the 
analysis using different exclusion masks, showing that a particular region of the sky 
near the Galactic plane shows a higher variance than 95.58% of the simulations whereas 
the rest of the sky has a lower variance than 99.96% of the simulations. The relative 
difference in variance between both regions is bigger than in 99.64% of the simulations. 
This anisotropic distribution of power seems to be causing the anomaly since the model 
assumes isotropy. Furthermore, this region has a clear frequency dependence between 
41GHz and 61GHz or 94GHz suggesting that Galactic foreground residuals could be 
responsible for the anomaly. Moreover, removing the quadrupole and the octopole 
from data and simulations the anomaly disappears. The variance anomaly and the 
previously reported quadrupole and octopole alignment seem therefore to be related 
and could have a common origin. We discuss different possible causes and Galactic 
foreground residuals seem to be the most likely one. These residuals would affect the 
estimation of the angular power spectrum from the WMAP data, which is used to 
generate Gaussian simulations, giving rise to an inconsistency between the estimated 
and expected CMB variance. If the presence of residuals is confirmed, the estimation 
of the cosmological parameters could be affected. 
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1 INTRODUCTION 

The Cosmic Microwave Background (CMB) is the most an- 
cient image of the Universe. This image can reveal valu- 
able information on the origin and evolution of the Uni- 
verse. Many observations in different cosmological fields tend 
to support the concordance, standard inflationary model. 
This model predicts the temperature anisotropics of the 
CMB to represent a Gaussian and isotropic random field on 
the sphere, whereas non-standard models, such as topologi- 
cal defect models, predict departures from isotropy and/or 
Gaussianity. 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite, provided the first all-sky, high resolution map 
(Bennett et al. 20031 ) of the CMB anisotropics. In the past 
years, many authors claimed to have found anomalies in 
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the WMAP data. Some of them could be affected by a 
posteriori statistics but the list of anomalies is too long 
to be ignored without worry. Some of these anomalies are, 
asymmetry between ecliptic hemispheres (JEriksen et al 
'2004'; 'Hansen et al .11 20041: lEriksen et a l."2004': 'Hansen et al 



2004; Eriksen et_al.l l2QQ5l: IPonoghu e & Donoghue 2 0oi: 
Land & Maguei^o 2005; Be rnui et al.1 I2Q07I: lEriksen et al 
2007 : Bernui 2009: Hoftuft et al.1 l2009l : IPietrobon et al 
bold : IVielva fc Sa^d |201QI : ICavonI |201Q|): anJ^- 
lous quadrupole- o ctupo le alignment rBielew icz et al.l 



| 2004l: ICopi et al.1 [20o3: IdeOliveira- Co sta et a l.1 l2004l: 
Schw arz et al.1 l2004l: iBielewicz et al.1 l2005l: lAbramo et al 



2006: ICopi et al.1 l2006l. l2007l : iGruppuso k Buriganal [20091 : 
Frommert fc EnfilinI |201QI 1 ; unexpected alignment of CMB 



structures towards the Ecliptic poles (|Wiaux et al.l l2006l : 
Vielva et al. 2007); a prominent cold s pot in the southern 



Gala c tic hemisphere (IVielva et al . 2004; Mukhe riee fc Wang 
l2004l : ICruz et al.1 l2005l . l2006l . l2007a .b.. ,2008,; .McEwen et al 
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20051 ; ICavon et al.1 [JOOsI: iRath et al.1 l2007l: IWiaux et aD 



20081: iGurzadvanfc Kocharvanllioosi: IPietrobon et al.ll2008l : 



Gurzadvan et al."2009': 'Rossmanit h et al.l l200^ 

These anomalies arise after using very sophisticated sta- 
tistical methods. Howe ver, a very simple analysis already 
reveals some anomalies. iMonteserm et al.l (|2008i ) estimated 
the CMB variance of the 3-year WMAP data and com- 
pared it to the values obtained from Gaussian simulations. 
The p- value of the variance was 1.2%. They confirmed this 
result by comparin g with an analytic formula derived in 
ICavon et al.l (|l99lh . An analysis of the two ecliptic hemi- 
spheres revealed that the anomaly was mainly due to the low 
variance of the Northern ecliptic h emisphere. Si milar results 
were found in iLarson k, Wandeltl (^2004 ) and .Avaita et al.l 
(12003), analysing spot abundances. 

Here we propose a new method to estimate the CMB 
variance from data with anisotropic noise. In section 2 we re- 
view the WMAP data and exclusion mask s we use. We com- 
pare t he new method to the one used in iMonteserin et al.l 
(|2008l l in section 3 and apply it to the 5-year WMAP data 
in section 4. Finally we present the conclusions of this paper 
in section 5. 



2 WMAP DATA AND EXCLUSION MASKS 



The WMAP team (|Hinshaw et al.1 l2009l ) provides 
foreground-reduced maps for the Q (41 GHz), V (61 
GHz), and W (94 GHz) bands. We use the WMAP 
foreground cleaned VW combined map: 





Figure 1. Mollweide projections of four exclusion masks used 
in the analysis. From top to bottom and left to right we have 
the kq85, kq75, gc30 and gc40 masks. The excluded regions are 
plotted in black. 



In our Galactic foreground analysis we build the gp33 
and gplO masks which admit only some pixels near the 
Galactic plane. The description of these masks is given in 
section 4. 

The simulations are generated from the best fit angu- 
lar power spectrum estimated from the data by the WMAP 
team. Note that the kg 85 mask is used in the angular power 
spectrum estimation (jNolta et al.ll2009l ) . Isotropy and Gaus- 
sianity is assumed and the beam and noise properties are 
taken into account. 

Note that the kpO and k p2 masks given in the 3 -year 
WMAP data release used in IMonteserm et al.1 (|2008h. are 
analogous to the kq75 and kq85 masks but slightly different. 



Te(x) = yT,(x) i/;,(x). 



(1) 



where x = iO^cj)) gives the position in the sky and the in- 
dex i stands for the 6 differencing assemblies of the V and W 
bands, namely VI, V2 for the V band and Wl, W2, W3 and 
W4 for the W band. These are the data recommended for 
cosmological analyses, since their Galactic foreground con- 
tamination should be low. The noise weight Wi{yi) is defined 
as: 

Wi{yi) 



Wi(yi) 



E^=l^4x) 



Wi{yi) 






(2) 



where aoi is the noise dispersion per observation for each dif- 
ferencing assembly and Ni (x) is the number of observations 
made by the receiver i at the position in the sky x. 

Although the data are provided at a HEALPi^j 
([Gorski et al.l 1 20051 ) resolution of Uside = 512, we degrade 
the data to Uside — 256 since the smallest scales are dom- 
inated by noise. In addition, in order to avoid the strong 
contamination present at the Galactic plane and the emis- 
sion coming from bright extragalactic point sources, we will 
apply different exclusion masks. 

The WMAP team provides the kq85 and kq75 masks 
which leave approximately 85% and 75% of the sky un- 
masked. We construct two more masks by applying an ad- 
ditional ±30 and ±40 degrees Galactic latitude cut to the 
kq75 mask. Hereafter we will call these masks gc30 and gc40 
respectively. These four masks are shown in Figure [1] 



3 CMB VARIANCE ESTIMATION 

Assuming that the foreground contamination has been to- 
tally removed, our data, T(x), consist of two components, 
namely CMB, Tcmb(x), and instrumental noise, Tnoise(x). 



T (x) = TcMB (x) ± Tnoise (x) . 



(3) 



The standard cosmological model predicts Tcmb(x) to 
represent an isotropic Gaussian random field. Hence, the 
CMB variance, Gq^ is constant for any position of the sky, 
X, according to that model. Therefore 



Tcmb(x)- Ar(0,cro) 



(4) 



On the contrary, the noise is Gaussian, uncorrelated, 
but anisotropic. Hence the noise variance depends on the 
position on the sky, (j'^oisei^)^ but it can be calculated since 
the number of observations per pixel is known 

Tnoise (X)- AT (0,aLse(x)). (5) 

Note that the correlated WMAP noise w a s sho wn to 
be negligible in section 4.2 of IMonteserm et al.1 (|2008l ). The 
aim of the present paper is to estimate cfq^ comparing the 
obtained value to that of Gaussian simulations and checking 
the isotropy of the CMB. We use a new estimator but let us 
review first the estimator used in iMonteserin et al.l ([20 



3.1 The minimum Kolmogorov-Smirnov distance 
estimator 



http://healpix.jpl.nasa.gov 



The normalised CMB temperature is given by 
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IX (x) 



Te(> 



\/^§ + ^noise (X) 



(6) 



According to the description given above, tx(x) has 
normal distribution: 



Kx)~iV(0,l). 



(7) 



In order to estimate the CMB variance, we replace Gq 
by So, which we allow to vary in a sufficiently large range. 
In this manner we can find the value which minimizes the 
Kolmogorov-Smirnov distance to a normalised Gaussian. 
This value is our estimate for (Tq. Hence we can write the 
minimum Kolmogorov-Smirnov estimator (hereafter MKS 
estimator) as: 



ctq = So : min 



; {Fu {si) 



$ 



(8) 



where $ is the Cumulative distribution function of a 
standard normal and Fu (sq) is the empirical distribution 
function for the normalised temperature. 



x10"^ 



Unit variance estimator 




true a^ (mK^^ 



3.2 The unit variance estimator 

Here we propose an alternative way of estimating Gq. We 
replace cfq by a set of values Sq varying in a sufficiently 
large range and find the value where the variance of the 
normalised map, var iu (sq)) is closer to unity. This value 
of Sq is the estimated Gq: 



Sq : mm var 



(«(sg)) 



(9) 



Let us check how this unit variance estimator (hereafter 
UV estimator) performs compared to the MKS estimator. 



3.3 Comparing both estimators 

Two estimators can be compared in terms of their mean 
squared error, which can be written as the sum of the vari- 
ance and the squared bias. For an estimator T trying to 
estimate the true value we have that the mean squared 
error reads: 



m.s.e(T) = E [(T - Of] = var(T) + briO), 
where briO) is the bias: 
briO) =E[T] -0. 



(10) 



(11) 



Since we only have one sky we estimate the mean 
squared error, variance and bias using simulations. Each 
simulation has a different Gq , which can be calculated before 
adding the noise. The estimation of this quantity, ctq, is then 
obtained applying the estimator to the simulation including 
the noise. Hence we take T = ctq — ctq and ^ = 0. Note that 
T corresponds to the distance of each dot to the red line 
in Figure [2] We apply the MKS and UV estimators to 10^ 
simulations using the kq75 mask. The results are shown in 
Table [1] The UV estimator performs better due to its lower 
bias and variance. 

Since we have found that th e UV estimator p e rform s 
better than the MKS one used in iMonteserin et ahl (|2008[ ), 
let us repeat the variance analysis with the new tool. 
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Figure 2. Estimated versus true dispersion using the UV esti- 
mator (top) and the MKS estimator (bottom) for 10^ Gaussian 
simulations. The kq75 mask is applied in both cases. 



4 THE ANALYSIS 

Applying the UV estimator to the WMAP foreground 
cleaned VW combined map we find ctq = 7.18 x 10~^ mK^ 
which compared to 10^ simulations gives a p- value of 
0.31%. The kq75 exclusion mask was used in this analy- 
sis. This p- value is much lower than the 1.2% one found in 
iMonteserm et al.l |2008) using the 3-year QVW combined 
map, masked with the 3-year kpO mask. 

The significance of the anomaly is high. No a poste- 
riori choices have been made since the WMAP foreground 
cleaned VW combined map and the kq75 mask are recom- 
mended by the WMAP team for such analyses. 

However the cause of the anomaly is unknown. In the 
following subsections we perform a number of follow-up tests 
to further investigate its possible origin, but these a poste- 
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Table 1. Estimated mean squared error, variance and bias of T = (Jq—cFq for the Unit Variance (UV) and Minimum Kolmogorov-Smirnov 
(MKS) estimators. 10^^ simulations have been used, applying the kq75 mask. 



Estimator T = cTq - 


--1 


Mean Squared Error (mK^) 


Bias {mK'^) 


Variance of the estimator (mK^) 


UV 




8.72 X 10-11 


3.43 X 10-'^ 


8.71 X 10-11 


MKS 




1.23 X 10-8 


-3.19 X 10-5 


1.12 X 10-8 



riori tests are not valid to establish the significance of the 
anomaly. To make this clear, we will not consider the lower 
tail probabilities of the follow-up tests as p- values. 

Let us start applying different exclusion masks to anal- 
yse the possible cause of the anomaly. 



4.1 Applying different exclusion masks 

In Table [2] we list the results for the kq85, kpO, kq75, gc30 
and gc40 masks in increasing order of masked pixels (see sec- 
tion 2 for the definition of the masks). The corresponding 
histograms are shown in Figu re [3l Th e three year kpO ex- 
clusion mask used in Montes erin et al.l (2008) excluded less 
pixels than the kq75 mask used here and that does change 
in some extent the significance as can be seen in Table [21 
Note that the lower tail probability decreases as the mask 
grows reaching its minimum at the gc30 mask and increases 
slightly for the gc40 mask. 

The estimation of the variance and hence the variation 
of the lower tail probability is affected by the efficiency of the 
estimator. If we assume that the low variance was intrinsic 
to the all-sky CMB data, we would expect the lower tail 
probability to have its lowest value for the smallest mask (i.e. 
kq85) since the variance of the estimator is higher when less 
pixels are taken into account. Applying the UV estimator 
to 10^ simulations using the kq75 and gc40 masks, we find 
that the m.s.e of the estimator is 36% higher for the gc40 
mask respect to the kq75 one. The estimation is worse since 
less pixels are used. This has to be the reason for the higher 
lower tail probability for the gc40 mask compared to the 
gc30 one. Note that the histograms in Figure [3] are broader 
for bigger masks due to this effect. 

However increasing the admitted fraction of the sky 
from the gc30 to the kq85 mask the lower tail probabil- 
ity grows instead of decreasing, hence this can not be due 
to estimator efficiency. This trend rather suggests a higher 
variance in the data near the Galactic plane, which could 
perhaps be due to Galactic foregrounds. Before further in- 
vestigating this result, let us analyse the ecliptic he mispheres 
separately as was done in iMonteserin et al.l ((20081). 



4.2 Hemispherical analysis 

iMonteserin et al.l (|2008l l found that the significance of the 
low variance was higher in the northern ecliptic hemisphere. 
Therefore we also analyse both ecliptic hemispheres sep- 
arately, finding the same result but again with increased 
significance. Applying the kq85 mask, the northern ecliptic 
hemisphere shows a lower tail probability of 0.01% while 
the southern one is compatible with the simulations (lower 
tail probability 36.23%). A similar result is found analysing 
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Figure 3. Histograms of the variance of 10^ simulations, using 
4 different masks, namely kq85, kq75, gc30 and gc40. The red 
vertical line represents the value obtained for the data. 



Galactic, instead of ecliptic hemispheres. We find a lower 
tail probability of 0.03% for the Galactic north and 36.82% 
for the Galactic south (see Table [3]). 



4.3 Localising the anomalous region 

The most significant lower tail probability in Table [2] is ob- 
tained for the gc30 case. Therefore the pixels excluded by 
the gc30 mask seem to have a high variance. In order to 
study these pixels we build a complementary mask which 
admits only the pixels between the kq85 mask and the gc30 
one (i.e. all the pixels not excluded by the kq85 mask but 
excluded by the gc30 one, see Figure |4}. We call this mask 
gp33 since the admitted fraction of the sky is about 33% 
and is located near the Galactic plane. The first two lines 
in Table [3] show the result obtained for the gp33 mask com- 
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Table 2. Estimated CMB Variance from the data (third column) using 5 different masks, namely kq85, kpO, kq75, gc30 and gc40. The 
sky fraction admitted by each mask at Ugide = 256 is shown in the second column, and the lower tail probability of the data calculated 
from 10^ simulations is listed in the right column. 



Mask 


Sky Fraction Admitted 


(%) 


Estimated Variance 


(mK^) 


lower tail probability 


kq85 


81 






7.33 X 10-3 




1.57% 


kpO 


76 






7.26 X 10-2 




0.70% 


kq75 


71 






7.18 X 10-3 




0.31% 


gc30 


48 






7.01 X 10-3 




0.23% 


gc40 


34 






6.94 X 10-3 




0.41% 




WMAP low variance region 



Figure 4. Mollweide projections of the eight exclusion masks 
used in this section. From top to bottom and left to right we 
have the gc30, gp33. Galactic north (galN), Galactic south (galS), 
ecliptic north (eclN), ecliptic south (eclS), kq80 and gplO masks. 
The excluded regions are plotted in black. 




• 




WMAP high variance region 




pared to the gc30 one. The gp33 result is compatible with 
the simulations and completely different from those found in 
Table [2] In the previous subsection we analysed separately 
Galactic and ecliptic hemispheres. The results were similar 
to those obtained when studying the complementary regions 
allowed by the gp33 and gc30 masks as we show in rows 3 
to 6 in Table El 

Putting these three results together, it is clear that the 
intersection of the three high- variance regions (Galactic and 
ecliptic south, and gp33) could have an even higher variance. 
This intersection, i.e. the region of the sky admitted by all 
three masks, represents about 10% of the sky. In order to 
study this region we build the gplO mask admitting this 
small fraction of the sky and its complementary mask which 
we call kqSO. In Figure [5] we show the Mollweide projection 
of these two masks applied to the WMAP VW combined 
map. The estimated variances for these masks are listed in 
the last two lines of Table [3] The gplO region shows an 



Figure 5. Mollweide projection of the kq80 (top) and gplO (bot- 
tom) masks applied to the WMAP combined VW map. The ex- 
cluded regions are shown in grey. 



anomalously high variance, whereas the kqSO region presents 
an anomalously low variance. The histograms corresponding 
to Table [3] are plotted in Figure [6] 

The relative difference of estimated variances between 
the WMAP gplO and kqSO regions, is bigger than in 99.64% 
of the simulations. Therefore, the different variance between 
the region near the Galactic plane and the rest of the sky 
seems to be violating the cosmological principle of isotropy 
and could be the cause of the anomalous variance. Isotropy 
is one of the key assumptions mad e to estimate the angular 
power spectrum (jNolta et al.ll2009r ) and to generate the sim- 
ulations. Let us recall the process we follow in this analysis. 
We use the best fit d estimated from the data to generate 
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Table 3. Variance of the data (third column) using 4 different pairs of complementary masks, namely Galactic north (galN) vs. Galactic 
south (galS), ecliptic north (eclN) vs. ecliptic south (eclS), gc30 vs gp33 and kq80 vs gplO (see text for detailed description of these 
masks). The sky fraction admitted by each mask at Uside = 256 is shown in the second column, and the lower tail probability of the 
data calculated from 10^ simulations is listed in the right column. 



Mask 


Sky Fraction Admitted (%) 


Estimated Variance 


{ruK^) 


lower tail probability 


gc30 
gp33 


48 
33 


7.01 X 10-3 
7.79 X 10-3 




0.23% 
38.07% 


galN 
galS 


40 
41 


6.84 X 10-3 
7.81 X 10-3 




0.03% 

36.82% 


eclN 
eclS 


41 
40 


6.86 X 10-3 
7.81 X 10-3 




0.01% 
36.23% 


kq80 
gplO 


72 
9 


7.04 X 10-3 
9.39 X 10-3 




0.03% 

95.58% 
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Figure 6. Histograms of the variance of 10^ simulations, using 
4 pairs of complementary masks. Each pair is plotted in a dif- 
ferent row. From top to bottom we show gc30 versus gp33 mask. 
Galactic north versus Galactic south, ecliptic north versus ecliptic 
south and kq80 versus gplO mask. The red vertical lines represent 
the values obtained for the data. 



Gaussian and isotropic simulations, and then we compare 
the variance of the simulations again to the data. 

Since the kci 85 mask is used for estimating the Ci 
(jNolta et alj l2009l l and there is an excess of power in the 
gplO region, the isotropic simulations will have an excess of 



variance in the complementary kqSO region when compared 
to the data. 

Hence the anomaly could not be due to an isotropic 
low variance of the data over the whole sky, but rather to 
an anisotropic power distribution. 

The high variance in the Galactic plane region is a clear 
hint for Galactic foreground contamination, since adding 
any two independent emissions ends up in an increase of the 
variance. However we would need some other support for the 
Galactic foreground hypothesis, since the anisotropic vari- 
ance distribution could be casual. Contaminating Galactic 
foregrounds such as synchrotron, free-free and Galactic dust 
emission show a strong frequency dependence. Therefore, 
a frequency dependence of the variance in the gplO region 
would be proving the foreground hypothesis right. 



4.4 Individual frequency bands 

Here we perform the analysis on the Q (41 GHz), V (61 GHz) 
and W (94 GHz) frequency bands separately. We compare 
the estimated variance to lO"^ simulations for each band. 

We find a slight frequency dependence of the lower tail 
probability when using the gplO mask but the significance 
remains almost unchanged using the kqSO one. The lower tail 
probabilities for the Q, V and W bands are 0.02%, 0.01% 
and 0.01% respectively, using the kqSO mask; and 96.7%, 
95.4% and 95.7% applying the gplO mask. 

In order to check whether the frequency dependence of 
the variance is significant, we calculate the difference of the 
estimated variance between different bands, and compare 
their histograms to the data in Figure [71 

The gplO variance differences including the Q band 
data, show a very high frequency dependende. This result 
indicates the presence of some foreground residuals near the 
Galactic plane in the Q band. The V and W bands which 
are the cosmological bands, do not show any anomalous fre- 
quency dependence, but the Q band result and the anoma- 
lously high variance make this region highly suspicious of be- 
ing contaminated. The combination of dust and synchrotron 
or free-free emission could perhaps give a flat frequency de- 
pendence at the V and W frequencies. Moreover an artifact 
of the foreground cleaning process could be affecting the 
pixels in this region. 
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CMB Variance Q - CMB Variance V, kq80 mask 
300 r 
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Figure 7. Histograms of the differences of estimated variances 
of the Q and V bands (top row), Q and W bands (middle row), 
V and W bands (bottom row). The kqSO mask has been used 
for panels of the left column and the gplO one for the right col- 
umn. 10^ simulations have been used for each band. The data are 
represented by a red vertical line. 
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Figure 8. Histogram of the variance removing the quadrupole 
and octopole from data and simulations. Four different masks are 
used. From top to bottom we apply the kq85, kq75, gc30 and gc40 
masks. 



The excess of variance near the Galactic plane could 
be related to other, previously reported anomalies. For in- 
stance, many authors (see introduction) reported an anoma- 
lous quadrupole-octopole alignment in a direction near the 
Galactic plane. 



out foreground residuals as the cause of the quadrupole- 
octopole alignment. The results in this section show us that 
if there are Galactic foregrounds in the gplO region, they 
mainly affect the quadrupole and octopole. 



4.5 Subtracting Quadrupole and Octopole 

Let us study a possible relation between the low variance 
and the anomalous quadrupole-octopole alignment. In or- 
der to check the possible influence of these multipoles in 
our results, we remove the quadrupole and octopole from 
data and simulations, and repeat the analysis. The results 
plotted in Figure [8] and Figure [9] show that the data are 
now fully compatible with the simulations. The link between 
both anomalies is clear and they could have a common cause. 
Note that Figure [8] still reveals a slight dependence of the 
variance with the mask, i.e. the less pixels in the Galactic 
region are admitted, the lower is the variance of the data as 
compared to simulations. This could hint again to the pres- 
ence of residual Galactic foreground contamination, but the 
va riation is not statistica lly significant. Some auth ors, such 
as ISlosar fc SeliakI (|2004l ) or lAbramo et all (|200 9V pointed 



4.6 WMAP 7 year data 

At the same time this paper was written, the WMAP team 
relea sed a new version of the WMAP data (jJarosik et alj 
120101 ) . The data and masks given in the new release are very 
similar to the 5-year ones. We estimated the CMB variance 
for the 7- year data using the 7- year kq75 and kq85 masks ob- 
taining 7.17 X 10~^mi^^ and 7.27 x 10~^mK^ respectively, 
which differ less than 0.1% and 0.9% from the 5-year re- 
sults listed in Table [2l Similar differences are found using 
the other masks. The lower instrumental noise of the 7- year 
data has a negligible effect in our variance estimation. The 
estimated variance for the 5-year data using the 7- year kq85 
mask is 7.26 x 10~^raK^ which is a variation of 0.1% with 
respect to the 7- year data result. Hence our results and con- 
clusions remain unaffected by these small changes. 
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Figure 9. Histogram of the variance removing the quadrupole 
and octopole from data and simulations. The kqSO (top) and gplO 
(bottom) masks are used. 



5 CONCLUSIONS 

In this paper we compare the variance of the WMAP data 
to Gaussian simulations generated from the best fit an- 
gular power spectrum estimated from the data. Isotropy 
and Gaussi anity are as sumed to generate the simulations. 
iMonteserin e t al. (2008) already performed this comparison 
for the 3-year WMAP data, using the MKS variance esti- 
mator. They found that the variance of the data was lower 
than 98.8% of the simulations. 

Here we apply the UV estimator to the 5-year WMAP 
data. The UV estimator performs better than the MKS one 
due to a lower estimator variance and bias. We find again a 
low CMB variance with increased significance. The p- value 
for the kq75 mask is 0.31%. Performing the same analysis ap- 
plying different exclusion masks we find a strong dependence 
of the variance with the used exclusion mask. The region 
near the Galactic plane shows a higher variance and com- 
patibility with simulations, while the complementary region 
near the Galactic poles shows an anomalously low variance. 
We further localise the high- variance region, finding that 
the excess in variance is caused by about 10% of the sky. 
This region, which we call the gplO region, is located near 
the Galactic plane and shows an anomalously high variance 
(lower tail probability = 95.58%), whereas its complemen- 



tary region (kq80 region) has an anomalously low variance 
(lower tail probability = 0.03%). 

The relative difference of estimated variances between 
the WMAP gplO and kq80 regions, is bigger than in 99.64% 
of the simulations. Hence, this anisotropic distribution of 
the variance could be the cause of the observed anomaly. 
The angular power spectrum has been estimated applying 
the kq85 mask, admitting the high- variance region near the 
Galactic plane. Since isotropy is assumed, this high power 
is distributed over the whole sky when generating simula- 
tions. Comparing this high-variance simulations to the low 
variance region of the sky, they turn out to be incompatible. 

Studying the relation of the low variance anomaly to 
the quadrupole-octopole alignment, we find that removing 
the quadrupole and octopole from data and simulations, the 
variance anomaly disappears. This link hints for a common 
cause of both anomalies. 

There are several possible causes for both anomalies, 
namely unknown systematic effects, fortuitous alignment of 
quadrupole and octopole. Galactic foregrounds and cosmic 
defects. Any unknown systematic effect could be causing 
the anomaly, although there is no evidence for this hypoth- 
esis yet. A fortuitous alignment of qu adrupole and octopol e 
has been proven to be very unlikely ([Bennett et al.ll201Ql ). 
A high variance in the region near the Galactic plane is a 
hint for Galactic foreground contamination. Furthermore we 
find strong evidence for for eground contami n ation in the Q 
band high- variance region. ISlosar fc SeliakI (| 20041 ) already 
point ed out Galactic foregr ounds as the cause of the align- 
ment. iBennetLeLajJ (|20ld ) shows that masking some small 
regions in the sky the alignment disappears. In particular, 
masking some spots in the high- variance region, eliminates 
any significant alignment. This shows that the alignment is 
not robust and could be due to Galactic foregrounds in the 
high- variance region. Topological defects such as textures 
(jCruz et al.ll2007bl ) could also be considered as an alterna- 
tive explanation. 

In our opinion the most likely hypothesis given the fre- 
quency dependence with respect to the Q band and the lo- 
cation of the high- variance region are Galactic foreground 
residuals. If this hypothesis was proven right, the cosmolog- 
ical parameter estimation could be affected in some extent. 
However, this has to be carefully analysed in future work. 
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